Excitation of coherent phonons in the one-dimensional Bi(114) surface 
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We present time-resolved photoemission experiments from a peculiar bismuth surface, Bi(114). 
The strong one-dimensional character of this surface is reflected in the Fermi surface, which consists 
of spin-polarized straight lines. Our results show that the depletion of the surface state and the 
population of the bulk conduction band after the initial optical excitation persist for very long times. 
The disequilibrium within the hot electron gas along with strong electron-phonon coupling cause 
the excitation of phonon standing waves, which, in turn are reflected in coherent modulations of 
the electronic states. Beside the well-known Ai g bulk phonon mode at 2.76 THz the time-resolved 
photoelectron spectra reveal a second mode at 0.72 THz which can be attributed to a standing 
optical phonon mode along the atomic rows of the Bi(114) surface. 



Reduced dimensionality at solid surfaces results in 
unique characteristics of the electronic structure as com- 
pared to the bulk. In the extreme case of a recently 
discovered phase of solids, the so-called topological insu- 
lators, the topology of the bulk band structure demands 
the existence of topologically protected metallic states 
at the solid surface [TJ [2]. Bismuth, a semi-metal on 
the verge of being a topological insulator, also features 
metallic surface states on all surfaces studied so far, giv- 
ing a strong enhancement of the metallic density of states 
(DOS) at the surface [5HS]- The degeneracy of the surface 
states in momentum-space is lifted due to the strong spin- 
orbit coupling in bismuth and the broken inversion sym- 
metry at the surface; this effect is called the Rashba-effect 
[BJ. In particular, the strongly anisotropic (114)-surface 
was found to support a spin-split quasi-one-dimensional 
(ID) metallic surface state [7J. The momenta \lf of the 
states at the Fermi energy Ep form straight lines close 
to the T-points of the surface reciprocal lattice and per- 
pendicular to the direction [110] of the atomic rows. In 
such a topology, any low-energy perturbation providing 
the momentum required to connect two parallel sections 
of the Fermi surface, generates a strong electronic re- 
sponse. The general response function can be calculated 
using the Linhard function. It becomes singular at mo- 
menta connecting long parallel sections of the Fermi sur- 
face. This leads to instabilities of the system, which may 
appear as charge- or spin-density waves, for instance [5]. 
Previously, such a charge-density wave (CDW) ground 
state was proposed for the Bi(lll)-surface [5]. On that 
surface the hexagonal shape of the Fermi contour of the 
electron pocket at T provides straight, parallel sections at 
opposite momenta ±kj?. A combined photoemission and 
tunneling microscopy study, however, provided no evi- 
dence for a CDW [10]. Kim and co-workers argued that 
due to spin-orbit interaction, the states at opposite mo- 



menta carry opposite spin, and any scattering between 
the states at these momenta requires spin-flip. 

Thus, the situation on Bi(114) is particularly inter- 
esting: The Fermi surface topology renders the metallic 
surface states instable against charge-density wave tran- 
sitions, because the complete Fermi surface of the surface 
states can be nested by one momentum vector. Moreover, 
the distorted rhombohedral lattice structure of bulk bis- 
muth makes it sensitive to displacive excitations of co- 
herent optical phonons (DECP) [TlTfTS] . Such coherent 
oscillations feed strong electronic modulations in space 
and time giving evidence for a particularly strong in- 
terplay of electron states and phonons. On the other 
hand, due to spin-orbit interaction, the spins of states 
at the momenta +kp and — have anti-parallel orien- 
tation [7], which inhibits scattering between these states 
and, thereby, the stabilization of a standing wave. The 
Bi(114) surface therefore offers the unique possibility to 
study the competing effects of spin-momentum locking, 
reduced dimensionality, electron-phonon coupling, and 
Fermi surface topology. 

The Bi(114) crystal was cleaned in situ by cycles of 
argon sputtering and annealing at 300 K [7J . The experi- 
mental pump-probe setup was described in detail else- 
where [16 . Briefly, we excite the sample with a lin- 
early p— polarized infrared pump pulse (hv\ — 1.55 eV). 
The duration of the pump pulse At pump could be tuned 
by means of a grating compressor. In the measure- 
ments shown here, the absorbed pump fluence was set to 
0.57 mJ/cm 2 , corresponding to an excitation density of 
n ss 0.34% of the valence electron density in bulk bismuth 
[T4"] . The electronic structure was probed by two-photon 
photoemission (2x hvi = 2x3.1 eV, At pro b e = 50 fs), and 
the photoelectrons were detected using a hemispherical 
electron analyzer |17j . The combined energy resolution 
of light and spectrometer was about 50 meV, the angu- 
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FIG. 1: (Colour online) Sketch of the atomic structure of (lx2)-reconstructed Bi(114); the atoms forming the rows are 
highlighted; bottom: Brillouin zones for surface and bulk in the plane of the binary axis and the surface normal; the projection 
of the bulk X-point along the binary axis corresponds to the center of the second surface BZ Tio. Note that the green shaded 
(110) plane is tilted out of the image plane. Panel (a)-(e): Photoemission data taken at the Fermi level with hv = 21.2 eV 
(a and b) and 2PPE with 2 x 3.1 eV (c and d): (a) and (c) Fermi surfaces; the surface BZ is indicated by a black rectangle. 
Black represents high intensity, (b) and (d) Two selected spectra taken for momenta indicated by the blue arrows in (c); the 
black open diamonds in (b) and (d) show the difference between the two spectra, highlighting the surface state at P. Note the 
different energy and momentum scales in the 2PPE data, (e) Band dispersion recorded using hp = 21.2 eV (spin-integrated). 
The splitting at Ep is clearly visible. The dashed lines serve as guides to the eye, the color code mimics the orthogonal spin 
polarisations. 



lar resolution was set to ±1°. All measurements were 
performed at room temperature. 

The Bi(114) surface consists of straight atomic rows 
running along [110], which are separated by 28 A wide 
valleys, as shown in Fig.[l][7]. The one-dimensional spa- 
tial character is reflected in the surface electronic struc- 
ture. The Fermi surface maps in Figs.JIJa) and (c) reveal 
straight lines parallel to TY, which were attributed to a 
spin-split ID surface state [TJ. All other features derive 
from bulk states. They disperse according to the bulk 
BZ and possess a defined parity and symmetry only with 
respect to the remaining mirror plane yz. The surface 
state is localized perpendicular to the atomic chains and 
strongly dispersing along the chains (TX). Energy spec- 
tra taken for various momenta along TX reveal the sur- 
face state on top of A-shaped bulk bands, and with the 
band bottom at a binding energy of 100 meV, as shown 
in Figs. [TJb) , (d), and (e). The agreement between the 
conventional one-photon photoemission and low-energy 
2PPE data, shown in Figs. [TJa)-(b) and (c)-(d), respec- 
tively, ensures that the low-energy data contain no contri- 
bution from intermediate states in the unoccupied regime 
between Ep and the vacuum level [15] . 

In Fig. [2] photoemission data are shown as func- 
tion of time delay, taken close to normal emission 
(k x = hp = 0.037 A -1 . Note that only the surface 
state contributes to the intensity at the Fermi level. At 
zero pump-probe delay and throughout the entire Bril- 
louin zone, the bulk conduction band is populated with 



a transient hot electron population following absorption 
of the infrared pump pulse. A few possible momentum- 
conserving optical transitions along [114] are depicted in 
Fig.[2]ja). The electrons excited into the conduction band 
relax by scattering and accumulate at the band bottom, 
where they appear as a broad feature in the spectra, 
Fig. [2jb).|2D] The evolution of the hot electron distri- 
bution can be seen in the false color plots in Figs. J2^c) 
and (d). 

The thermalization and, eventually, the energy dissi- 
pation last over several picoseconds. In order to investi- 
gate the dynamics of the hot electron gas quantitatively, 
spectra and transients were fitted using Fermi-Dirac dis- 
tributions and rate equations, respectively. The results 
are plotted in Fig. [3] The electronic temperature T e i 
rises up to about 2000 K with a time constant of 260 fs. 
The maximum temperature and the time constant are 
in agreement with results of a previous x-ray diffraction 
study [H] and estimates based on the fluence used here 
|22) . respectively. The initial rise is followed by a slow 
cooling of T e i and energy dissipation from the electronic 
system to the lattice within about 5.9 ps. 

The cooling of the electronic system is accompanied 
by a corresponding shift of the band bottom on the same 
timescales, as shown in Fig. (3^b): Ignoring at present the 
first spike-like increase of the energy position, the aver- 
age peak position shifts by 60 meV towards lower energies 
within a few hundreds of femtoseconds and subsequently 
relaxes back on a picosecond timescale. Similar obser- 
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FIG. 2: (Color online) Time-resolved photoemission data: (a) Momentum conserving transitions from the 5th to the 6th 
valence band along the [114]-direction, which can contribute to the hot electron population observed close to F; the green lines 
denote the energy bands calculated in a tight-binding scheme [19]; electrons are excited by means of 1.55 eV photons (short red 
arrows) and accumulate in a minimum of the conduction band (grey shaded energy range), where they are probed by 2PPE 
(blue arrows), (b) Photoelectron spectra (symbols) taken close to F for various pump probe delays between -0.05 ps and 4 ps 
showing the pronounced signature of hot electrons above Ef ; The thick (thin) orange line denotes an exemplary fit using two 
Gaussians multiplied with a Fermi-Dirac distribution (bare Fermi-Dirac distribution) . (c) Plot of the photoemission intensity 
as function of energy and time delay, (d) Same data after subtraction of the hot Fermi-Dirac distribution for each individual 
delay in order to highlight the transient occupation of the conduction band. The black line follows the transient energy position 
of the bulk conduction band. 



vations were made very recently on Bi(lll): following 
absorption of an infrared pump pulse both, surface and 
bulk states shifted to lower energy [23] . The transient fol- 
lows the electronic temperature like in our case, and the 
authors showed that the shift in energy is caused by an 
electronic effect. Since the origin could not be elucidated 
so far, we conjecture that the shift is caused by a charge 
redistribution due to the excitation of hot electrons. 

The intensity cross correlation curves reveal two 
timescales Tj and Tu depending on the electron energy, 
as can be seen in Figs, ^c) and (d): The long decay 
time Tu for electrons from the conduction band bot- 
tom or below reflect thermally excited states for delays 
larger than 1.5 ps and, thereby, the evolution of T e i. 
The faster timescale for energies within the conduction 
band (about 0.6 eV and higher) is caused by comparably 
fast scattering within the conduction band, which domi- 
nates the decay of these states. The energy dependence 
corresponding to a power law with an exponent of —1.5 
depends on details of the transient electronic distribution 
function [24]. A detailed discussion is beyond the scope 
of this paper. 

The dynamics found here for Bi(114) strongly resemble 
those observed very recently in the topolocial insulator 
Bi 2 Se 3 [2S]: a fast intraband decay by electron-phonon 
scattering within the conduction band is followed by a 
much slower decay on a picosecond timescale for both the 
bulk conduction band and the metallic surface state. The 
second, long timescale is a consequence of the low scat- 
tering probability due to the low DOS at and above Ep in 
Bi on one hand, and due to the energy difference between 



conduction band and surface state, which is large com- 
pared to typical phonon frequencies. In contrast to the 
above mentioned study, however, we find strong evidence 
for simultaneous emission of many quanta of particular 
phonon modes from the observation of intensity and en- 
ergy modulations of the photoelectron peaks, which can 
readily be recognized in Fig. [5Jc) . 

In order to quantify these effects, the transient inten- 
sity changes were obtained by recording the intensity 
close to Ep in case of the surface state and of the conduc- 
tion band. The data were normalized to those at negative 
delay times. Two selected transients are displayed for two 
different pump pulse durations in Fig. [4ja) . The cosine- 
like modulations, shown in Fig. [4^b) , are assigned to the 
excitations of coherent phonon modes. Such modes are 
well known to modulate the bulk and surface electronic 
states, allowing them to be directly observed in photoe- 
mission spectra [2l)] - f2T)] : the excitation of many electronic 
transitions and strong coupling to specific phonon modes 
drives phonons anharmonically, which results in a crys- 
tal structure closer to the simple cubic lattice in the 
present case. The atoms then start to perform damped 
oscillations around the new equilibrium positions. This 
phenomenon is called displacive excitation of coherent 
phonons if the initial excitation occurs within a time pe- 
riod shorter than half a phonon cycle. As a consequence, 
the electronic spectral weight and energy position follow 
the changing atomic positions. Given that the electronic 
response is faster than the time period of the phonon 
mode, the phase of the coherent phonon is directly im- 
printed in the photoemission spectra. The amplitude of 
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FIG. 3: Analysis of hot electron dynamics: (a) Hot electron 
temperature T e i resulting from fits and (b) transient position 
of the conduction band as function of pump-probe delay, (c) 
Cross correlation curve obtained from the conduction band 
at 0.57 eV above Ef [Fig. [2jd)]; the solid line corresponds 
to a double-exponential fit with two decay constants n and 
m\ the dashed lines represent the two distinct contributions, 
(d) tu and r, as obtained from the fits as function of energy 
above Ef- The line represents a power law fit, which yields 
an exponent (E — Ep)~ lm6 '. 



the oscillations in binding energy (Fig. |3p) is of the or- 
der of 20 meV, which is in agreement with photoemis- 
sion data taken recently from Bi(lll) under comparable 
conditions and corresponds to an atomic displacement of 
about 1-2 pm (23] . 

The Fourier transform (FT) of the measured oscilla- 
tions, displayed as the top trace in Fig. [4jc) for a tempo- 
ral pump pulse width At pump of 160 fs, reveals two dom- 
inant frequencies: The first frequency of 2.76(2) THz is 
close to the softened longitudinal optical (LO) A lg mode, 
found at roughly 2.85 THz in numerous time-resolved 
experiments [TTHTS] at the T-point of the bulk BZ. The 
second frequency of 0.72(1) THz corresponds to an os- 
cillation period of 1.39(1) ps. The observation of this 
mode can be explained by the peculiar band structure of 
Bi(114) as will be discussed below. 

Assuming the electron-phonon coupling strength to 
scale with inverse phonon frequency, the FT may be 
compared to F{uj)duj/oj [3T], where F(ui)duj denotes the 
phonon DOS in the frequency interval [l),u + duS\. This 
function, which was calculated using the phonon DOS 
for two different electron excitation densities [30], is dis- 
played in Fig.gc). Two van Hove singularities dominate 
these spectra, one around 2.8 THz for the optical modes 
at r, and the second at about 0.7 THz corresponding to 
modes at the zone boundary along TKX of the bulk BZ. 
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FIG. 4: (Color online) Analysis of the modulations of photoe- 
mission intensities, (a) Transient photoemission intensity at 
the conduction band for a short pump pulse (top panel) and 
at Ef for a long pump pulse (bottom panel) as function of 
time delay. The solid lines denote rate equation fits, (b) As 
in (a) but after subtraction of the fitted rate equations: the 
oscillations of the intensity are clearly visible; the solid lines 
are cosine-functions at the dominant frequencies, (c) Fourier 
transform of the measured intensity cross correlations for the 
two different temporal pump pulse durations (solid symbols: 

^tvurr 



160 fs; open symbols: At purl 



280 fs). Bottom 



traces: calculated phonon DOS F(cj)du, taken from Ref. 1301 
and weighted by for the electronic ground state (blue) 
and for n = 1% excited valence electrons (red). 



We focus on the low- frequency mode at 0.72 THz. The 
excitation of this mode can either be the result of a strong 
decay channel of the high-frequency Ai g mode or else re- 
sult itself from a coherent electronic excitation. In order 
to discriminate between both excitation pathways, the 
pump pulse duration was increased to At — 280 fs by 
introducing linear chirp by means of a grating compres- 
sor. Since half the temporal period of the A\ g mode 
corresponds to about T^i g /2 s» 180 fs, the condition 
At < Tjiig/2 for coherent phonon excitation is no longer 
fulfilled. As a consequence, the A\ g mode is suppressed 
[see Fig. |4jc)]. The low-frequency mode at 0.72 THz, on 
the other hand, persists giving evidence for a coherent 
excitation of this mode via electronic transitions. 

The fact that this mode can very efficiently be directly 
excited in Bi(114) can be explained by the low dimen- 
sionality of the surface [5). The strongly enhanced sus- 
ceptibility of surface state electrons to excitations with 
momenta along the atomic chains increases the coupling 
significantly. Indeed, the strong phonon DOS around 
0.7 THz in Fig. gc) is due to van Hove singularities 
in three branches midway between two adjacent recip- 
rocal lattice points along TKX. This corresponds to a 
standing wave in real space along the binary axis [110], 
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i.e. along the atomic chains of the Bi(114)-surface. The 
character of these branches is acoustic in the long wave- 
length limit and changes to two transverse TO(X) and 
one longitudinal optical LO(X) branch at the border of 
the BZ[30 . Previously, a weak signal of about 1/6 of 
the Ai g intensity at 0.68 THz was reported from time- 
resolved reflectivity measurements on Bi films [32] and 
interpreted as transverse-acoustic mode at the X-point 
of the BZ in agreement with our interpretation of the 
Fourier spectra in the present work. Surprisingly and in 
contrast to a previous study on Bi(lll) [23], the surface 
states are strongly modulated. This indicates that we 
are indeed dealing here with a strong phonon mode in 
the chains of Bi(114). 

The wavelength of 4.54 A of these standing waves 
equals the atomic spacing in these chains, because the 
projection of the bulk X-point onto k x corresponds to 
the center of the second surface BZ as sketched in Fig. [T] 
This in turn means that the periodic lattice distortion 
must comprise at least the two topmost atomic rows of 
the chains (see Fig. [I]), and that the atoms in the ba- 
sis of the surface unit cell oscillate with respect to each 
other, which is equivalent to an optical phonon with in- 
finite wavelength. So far no evidence was found in our 
time-resolved experiments for phonon modes introduc- 
ing new spatial periodicities along the atomic rows. The 
spin-orbit interaction, which lifts the spin degeneracy at 
the surface, protects the metallic surface states against 
a phonon-driven charge-density wave transition, because 
efficient Fermi surface nesting is only allowed between 
states of the same spin helicity as shown convingly very 
recently for Pb wires on vicinal Si(557) [55] . 

In conclusion, the transient occupation of the bulk con- 
duction band of Bi leads to the excitation of at least two 
coherent optical phonon modes in the Bi(114) surface. 
Beside the A\ g bulk phonon a second mode at 0.72 THz 
is found to strongly modulate the surface electronic states 
on a picosecond time scale. This mode could be identi- 
fied as an electronically induced coherent displacive ex- 
citation of a standing wave along the atomic rows of the 
vicinal Bi(114) surface. 
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